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The Walecka many-body field theory is investigated in the context of quantum nonextensive 
statistical mechanics, characterized by a dimensionless parameter q. In this paper, we consider 
nuclear matter described statistically by a power-law distribution which generalizes the standard 
Fermi-Dirac distribution {q = 1). We show that the scalar and vector meson fields become more 
intense due to the nonextensive effects (g 7^ 1). From a numerical treatment, we also show that 
as the nonextensive parameter q increases, the nucleon effective mass diminishes and the equation 
of state becomes stiffer. Finally, the usual Maxwell construction seems not to be necessary for 
isotherms with temperatures in the range 14Mev < fcsT < 20MeV. 
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I. INTRODUCTION 

In the past few years, considerable attention has been 
paid to the so-called nonextensive statistics, both from 
theoretical and observational viewpoints. The nonexten- 
sive framework, first developed by Tsallis [l|, Q , seems to 
present a consistent theoretical tool to investigate com- 
plex systems in their nonequilibrium stationary states. 
Recently, several consequences (in different branches) of 
the Tsallis framework have been investigated in the lit- 
erature 0], which includes systems of interest in high 
energy physics^. In this regard, the very first applica- 
tion was addressed to a problem of solar neutrino [4|] . In 
this study, the authors applied the g-framework to derive 
a distribution function for the interior plasma, which is 
relevant in the calculation of the nuclear fusion reaction 
rates responsible for the neutrino flux emitted by Sun. 
Recently, within Fokker-Planck dynamics, the Einstein's 
relation between drag, diffusion, and the equilibrium dis- 
tribution for a spatially homogeneous system have been 
generalized in the context of Tsallis statistics, and such a 
model was applied to charm quark dynamics in a thermal 
quark^luon plasma for the case of coUisional equilibra- 
tion [5|. More recently, interpretations for central Au- 
Au collisions at RHIC energies in a Relativistic Diffusion 
Model (RDM) for multiparticle interactions based on the 
interplay of nonequilibrium and local equilibrium were 
studied also through extensive and nonextensive statis- 
tics 01 . 

In the theoretical treatment of the properties of nuclear 
matter, the relativistic phenomenological approach devel- 
oped by Walecka (the so-called quantum hadrodynam- 
ics [QHD-I]), represents one of the important approaches 
to the highly nonlinear behavior of strong interactions at 
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the hadronic energy scales. This model provides a ther- 
modynamically consistent theoretical framework for the 
description of bulk static properties of strong interacting 
many-body nuclear systems^. 



In this work we investigate the effects of the nonex- 
tensive statistical mehanics on the QHD-I theory by con- 
sidering a g-generalization of the Fermi-Dirac (FD) dis- 
tributions. In this regard, it is worth emphazing that 
earlier generalizations of quantum statistics for fcrmions 
and bosons have been developed in Ref. ,10,] and ap- 
plied in the context of relativistic nuclear equation of 
state (EoS) of Boguta-Bodmer [ll| by Drago et al p^ . 
However, in the present work, to study the nonexten- 
sive effects on the QHD-I theory we use the most recent 
fermion distribution obtained by Tweldeberhan, Plastino 
and Miller in Ref. ^ (TPM distribution). This devel- 
opment differs from the above cited in that it provides a 
new cut-off prescription based on the extremization of a 
thermodynamical functional. This paper is organized as 
follows. In Sec. II, we present the basic formalism of the 
mean field theory of QHD-I for the calculation of some 
important quantities of nuclear matter. A brief review of 
Tsallis statistics and TPM distribution function is made 
in Sec. III. In Sec. IV the convergence of the calculation 
is also considered and it is shown that the allowed values 
of q lie in the range 1 < g < 5/4. Our main results are 
discussed in Sec. V. We summarize the main conclusions 
in Sec. VI. 



^ As is well known, although providing a consistent theoretical 
framework, QHD-I presents some limitations. In this regard, a 
more complete theoretical treatment (for zero temperature) was 
discussed in Ref. whereas the thermo field dynamics in hot 
nuclear matter was considered in Ref. 0]. Here, however, we limit 
ourselves to the effects of nonextensivity on the QHD-I theory, 
as originally developed in Ref. 0]. 
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II. BASICS OF QHD-I 

As widely known, the Lagrangian density describing 
the nuclear matter reads [T] 



Additionally, we also use for the coupling constants the 
values of reference Q, namely'^, 



^ \ = 11.798 fm2 and f ] = 8.653 fm^ 



(9) 



which are fixed to give the bind energy -Ebind — —15.75 
MeV and kp = 1.42 fm^\ 



which represents nuclear matter composed by nucleons 
coupled to two mesons, namely, the cr and uj mesons (for 
details see reference Q). 

Applying standard techniques from field theory and 
the mean-field approach, we obtain the scalar density 



QS 



7n 



M* 



(27r)3 J E*{k) 
where M* is the effective mass 

M* = M - g„a = M 



[n{iy,T)+n{iy,T)]d^k , (2) 



rps 



(3) 



The baryon number density, the energy density and pres- 
sure are given, respectively, by 
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[n{i',T) - n{v, T)]dPk, 
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where 



and 
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n(v, T) = 



(M - M* f + 
n{v, T) + n{v,T)](fk 



,j3[E'{k)-u] I 



n{v,T) = 



(6) 



(7) 



(8) 



are the usual FD distributions for baryons and anti- 
baryons, with E*{k) = \/k^+M*^, (3 = l/ksT. The 
parameter v = fi — g^ujQ = M — [g^ IfnujY Qb is the effec- 
tive chemical potential, and 7a' is the multiplicity factor 
(7Ar = 2 for pure neutron matter and = 4 for nuclear 
matter) . 



III. NON-EXTENSIVE FRAMEWORK 

Over the past two decades, nonextensive statistical 
mechanics has successfully addressed a wide range of 
nonequilibrium phenomena in non-ergodic and other 
complex systems [1, 2]. As widely known, nonextensive 
statistical mechanics, as proposed by Tsallis [1], is a gen- 
eralization of Boltzmann-Gibbs (BG) statistical mechan- 
ics based on the functional 



w 

S, = -kBj2-' 



p! lllq Pi , Sq=l = 



W 

'ks'^Pi In- Pi 



(10) 



1=1 



Here, q is the nonextensive parameter and the q- 
logarithmic function above is defined as 



P 



1-9 _1 



1-g 



(11) 



In this g-framework, the additivity for two probabilis- 
tically independent subsystems A and B is generalized 
by the following pseudoadditivity 



S: 



AUB 



= + 5f + (1 - q)S^S' 



(12) 



where the cases q < I and q > I correspond to super- 
additivity and sub-additivity, respectively. For subsys- 
tems that have special probability correlations, extensiv- 
ity is not valid for Boltzmann-Gibbs entropy, but may 
occur for Sq with a particular value of the index q ^ 1. 
Such systems are sometimes referred to as nonextensive 



A. Quantum Statistics 

We recall the main aspects of the connections between 
the quantum statistics and Tsallis framework. Specifi- 
cally, we concentrate on the new cut-off prescription re- 
lated to the Tsallis' maximum entropy distributions (for 
details see Ref. [l^). The main result in this study is 



^ For the purpose of the present work, the values given in Eq. !9l 
suffices to investigate the effects of the nonextensivity in neu- 
tron and nuclear matter. Variations of the coupling constants, 
within the acceptable values given in current literature, do not 
qualitatively affect the conclusions. 
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FIG. 1: The effective nucleon mass and tfie vector meson field of pure neutron matter (7jv = 2) as function of temperature 
for different values of tlie parameter q. Panel (a): the self-consistent nucleon mass at vanishing baryon density. Panel (b): the 
vector meson field at nonzero baryon density corresponding to v = 540 MeV. 



that, for (7 > 1, a g-generalized quantum distributions for 
fermions and bosons are given by 



1 



e,(/3(e-A*))±l' 



(13) 



where reads 



eg{x) 



[1 + {q - l)x]i-^ if x>0 
[l + {l-q)x]^ if a;<0 



(14) 



and X — [3{e — /i). In the q ^ 1 limit, the standard 
FD distribution, n{iJ,,T), is recovered. As physically ex- 
pected, as T ^ 0, nq{^,T) n{^,T). This amounts 
to saying that for studies of the interior of neutron stars 
(where, in nuclear scale, T ~ 0) we do not expect any 
nonextensive signature. On the other hand, in heavy 
ions collision experiments or in the interior of protoneu- 
tron stars, with typical stellar temperatures of several 
tens of MeV (1 MeV= 1.1065 x W° K), nonextensive 
effects may appear. In order to study such effects, in the 
next Section we combine figs.©-® with the generalized 
FD distributions given by Eqs. and ([1] 



IV. NONEXTENSIVITY AND QHD-I 

Since the above function eq{x) is a deformed expo- 
nential, we must verify the mathematical convergence 
of the integrals in Eqs. ([3])- ([6]) by considering the q- 
distributions. Thus, let us write the eq{x) in a more 



compact form given as 

~eq{x)^[l + ^xf^ (15) 
where ^ = (7 — 1. Now, taking into account that 

~eq{x)^e'^x^'^. (16) 

in the limit x >> 1, we obtain the asymptotic behaviour 
for the integrals appearing in Eqs. ([3])- ([6]), i.e.. 



M* 
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M* d^k 



E*{k){eq[E*{k)-u]+l} 
d^k 



PB 
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-\E*{k)~v] + l 



E*{k) d^k 
[E*{k) + l 



fc2 d^k 



k' 



E*{k){eq[E*{k)^i^] + l]} ^1/? 



(17) 
(18) 
(19) 
(20) 



From Eqs. (fT7|l - ([20|) the general asymptotic behaviour 
can be summarized by fc^//c^/'>, for A^ = 2, 3, 4. In order 
to have k^ /k^^^ — > 0, when k — > 00, we find that 
1/^ > A, from which we obtain 



q < 



A+ 1 



A 



(for q> 1) 



(21) 



Note that, to simultaneously satisfy the convergence of 
all integrals in Eqs.©-®, we find that 1 < g < 5/4, 



4 




FIG. 2: Isotherms of nuclear matter (7jv = 4) equation of state at finite temperatures for different values of the parameter q. 



consistent with the limits of q in the range q G (0, 2] of 
Ref. dl. 

At high temperatures (T — > oo), the analytic solution 
to Eq. ((21) can be written as 



M 



(22) 



where £,i{q) is given by Eq. ((27)) . 

Several limiting cases of the EoS are of interest: 

1. The baryon distribution becomes a step function 
nq{v,Q) — 9{kp ~ |k|) in the limit T ^ for any 
value of . 

2. The system becomes degenerate in the limit ps ^ 
oo at any T. 

3. For ksT « M and pB 0, the equation of state 
of a classical nonrelativistic gas is obtained: 



e = [(2 - q)MpB + IpBksT] (23) 

7 — bq 2 

and 



P = 



7 -bq 



[{q-l)AIpB + PBkBT] . (24) 



Thus, the equation of state of a nonrelativistic gas 
of baryons can also be given by 

2{q-l)M + 2kBT 



^ ~ [2{2 - q)M + SkBT] ^ ■ ^^^^ 

Note that, as q 1, the limits of Ref. 0] are 
recovered. 



4. As r — !■ oo for any vaue of pB, an equation of state 
similar to that of a black body is obtained: 



In 



a{q)ikBT)\ p^e/3 



where, in Eqs.dUl) and 

z" dz 



eq{z) + 1 



(26) 



(27) 



When q ^ 1 we have that ^i{q ^ 1) — > 7r^/12 and 
£2(5 ^ 1) — > 77r''/120, recovering the limits of Ref. 



V. RESULTS 

In order to study the nonextensive effects on QHD-I 
theory, we have calculated numerically, for several values 
of temperatures and of the parameter q, the effective nu- 
cleon mass, the vector and scalar mesons fields for pure 
neutron matter {jpf = 2), as well as the EoS for symmet- 
ric nuclear matter {jn = 4). 

In Fig. [T]we show the nonextensive effects on the ef- 
fective mass and vector meson for a pure neutron mat- 
ter. In Panel (a), for = = 0, we note that the 
higher the parameter q, the smaller the effective mass 
M* (and, consequently, the higher the scalar field cr, 
since gacr = M — M*). Such an effect may be physically 
understood in that at a given temperature, the scalar 
density, as source for the scalar mesons, increases with 
the increasing of the nonextensive parameter q. Thus, 
the attraction of the nucleons, mediated by the scalar 
mesons, becomes stronger, reducing the effective mass. 
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FIG. 3: Panel (a): Isotherms of nuclear matter {"fM = 4) at T = 10 MeV for several values of the parameter q. Panel (b): 
Isotherms of nuclear matter corresponding to the critical temperatures and parameters Qc given in the data. The curves and 
the data are in the same order from top to bottom. Panel (c): the critical parameter as funtion of the critical temperature 
Tc corresponding to the data of Panel (b). 



The behaviour of vector meson field ^i^wq = {gui/Tn^Y qb 
is shown in Panel (b). 

Fig. [5] shows the nonextensive effect on the symetric 
nuclear matter EoS. The panels display the logp — loge 
plane for selected values of q. The results are ploted for 
arbitrarily chosen values of temperature, T =15 MeV, 
20 MeV, and 30 MeV. In reality, the motivation for this 
choice is of astrophysical interest in the study of pro- 
toneutron stars. Clearly, the nonextensive effect is man- 
ifested in the increasing of the pressure with the values 
of q becoming the EoS stiffer. 

Another interesting effect of the nonextensivity on nu- 
clear EoS of QHD-I concerns the phase transitions. From 
Panels (a) and (b) of Fig. [51 we see that the first order 
phase transition may be eliminated by the variation of 
the parameter q. This fact can be easily visualized in 
Fig. [3{a), where isotherms at T = 10 MeV are plot- 
ted for several values of q G [1.0, 5/4). Note that, for 
increasing values of g, the dip in the region of thermo- 
dynamical instability becomes smaller, vanishing at the 
turning point that defines the critical values of thermo- 
dynamical quantities (Tc, Pc, etc.). We also note that for 
our choice T = 10 MeV the upper value of q (near the 
5/4 limit discussed earlier) is not sufficient to eliminate 
the first order phase transition 16]. On the other hand, 
the phase transitions can be eliminated for values of T 
in the region 14MeV < T < 20MeV. In this interval, 
all temperatures can be made critical. This amounts to 
say that a (critical) parameter q^ can be determined in 
order to yield a turning point in the isotherm at a given 
temperature. This is illustrated in Fig. [3Kb) where sev- 



eral isotherms are displayed. In Fig. [HJc) values of q^ are 
plotted as function of temperature in the above range 
(14 MeV < T < 20 MeV). The smooth behavior of q^ 
allows a parametrization of the form 

qc = 1.7724(±0.003) - 0.0384(±0.0001)-— |— . (28) 



VI. FINAL REMARKS 

In this paper, we have investigated the nonextensive 
effects on the mean field theory of Walecka (QHD-I) 
7]. We have used, instead of the standard Fermi-Dirac 
nucleon and antinucleon distribution functions, the q- 
quantum distribution recently obtained by Teweldeber- 
han, Plastino and Miller in Ref. [l3|, which has a new 
cut-off prescription based on the extremization of a ther- 
modynamical functional. We emphasize that the nonex- 
tensive effects on nuclear and pure neutron matter, for a 
considerable range of temperature, is to make the equa- 
tion of state stiffer and to increase the intensity of the 
vector and scalar meson fields, with a consequent low- 
ering of the nucleon effective mass (for increasing val- 
ues of the nonextensive parameter q) . We believe that it 
may have interesting consequencies in astrophysical stud- 
ies, mainly in what concerns the calculation of masses of 
compact objects, such as protoneutron stars. 

As discussed in Sec. V, another interesting feature 
of nonextensivity is that, at temperatures in the range 
14 MeV < r < 20 MeV, phase transitions of first order 
can be avoided by a convenient variation of the param- 
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eter q, which allows the determination of a critical qc 
parameter at the turning point of an isotherm at a given 
T [Eq. Finally, it is worth mentioning that the es- 

timates for the nonextensive parameter from the gener- 
alized Walecka many-body field theory is consistent with 
the upper limit q < 2, obtained from several independent 
studies in the quantum [15] and in the non quantum limit 



[17| involving the Tsallis nonextensive framework. 
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